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Abstr act. The first part of the course is devoted to the study of solutions to the Laplace 
equation in Q \ K , where Q is a two-dimensional smooth domain and K is a compact 
one-dimensional subset of Q . The solutions are required to satisfy a homogeneous Neu- 
mann boundary condition on K and a nonhomogeneous Dirichlet condition on (part of) 
dfl . The main result is the continuous dependence of the solution on K , with respect 
to the Hausdorff metric, provided that the number of connected components of K re- 
mains bounded. Classical examples show that the result is no longer true without this 
hypothesis. 

Using this stability result, the second part of the course develops a rigorous mathemat- 
ical formulation of a variational quasi-static model of the slow growth of brittle fractures, 
recently introduced by Francfort and Marigo. Starting from a discrete-time formulation, 
a more satisfactory continuous-time formulation is obtained, with full justification of the 
convergence arguments. 
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1. NEUMANN PROBLEMS IN DOMAINS WITH CRACKS 

In these lectures il is a fixed bounded connected open subset of R 2 with a Lipschitz 
boundary dfl, and /C(f2) is the set of all compact subsets of fi. Given K £ /C(fi), we 
consider the differential equation Au = on £l\K with a homogeneous Neumann boundary 
condition on dK and on a part 9/vf2 of the boundary of f2, and with a non- homogeneous 
Dirichlet boundary condition on the rest of the boundary of VL\K . For simplicity we assume 
that 9/vf2 is a (possibly empty) relatively open subset of d£l, with a finite number of 
connected components, and we set <9dJ1 := <9f2\(5/v^ , which turns out to be a relatively 
open subset of dQ , with a finite number of connected components. For the applications 
we have in mind K will be a one-dimensional set, and will be considered as a crack in the 
domain f2, but we do not need this assumption in this lecture. 

If dK is not smooth, the variational formulation of this problems requires the Deny-Lions 
space L 1,2 (A), defined for every open set A C M 2 as the space of functions u £ Lf oc (A) 
such that the distributional gradient Vu belongs to L 2 (A;M. 2 ) (see [17]). The advantage of 
this space is that the set {V« : u £ L 1:2 (A)} is closed in L 2 (A;M. 2 ) even if dA is irregular 
(for the proof we refer, e.g., to [25, Section 1.1.13]). It is well known that, if A is bounded 
and has Lipschitz boundary, then L 2 (A;M. 2 ) coincides with the usual Sobolev space 7J 1 (A) 
(see, e.g., [25, Corollary to Lemma 1.1.11]). 
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(1.1) 



Given K £ JC(Q) and g £ i 1,2 (Sl\i4T) , we consider the following boundary value problem: 

Au = inO\^, 
du 



h/ on d{Q\K) n (K U d N Q) 



(1.2) 



u = g on 9c \ if . 
By a solution of (1.1) we mean a function u which satisfies the following conditions: 
u £ L^ 2 (n\K) , u = g on d D Q\K , 

/ VuVzdx = VzeL 1 ' 2 ^), z = on d D Sl\K , 
Jn\K 

where the equalities on dpi! are in the sense of traces. 

It is clear that problem (1.2) can be solved separately in each connected component 
of Q\K . By the Lax-Milgram lemma there exists a unique solution in those components 
whose boundary meets 8d^\K , while on the other components the solution is given by an 
arbitrary constant. Thus the solution is not unique, if there is a connected component whose 
boundary does not meet do^\K . Note, however, that Vw is always unique. Moreover, the 
map g i— > Vu is linear from L 1 ' 2 (Q\K) into L 2 (i}\K;R 2 ) and satisfies the estimate 

(1.3) / \Vu\ 2 dx< [ \Vg\ 2 dx. 

Jn\K Jn\K 

By standard arguments on the minimization of quadratic forms it is easy to see that u 
is a solution of problem (1.2) if and only if u is a solution of the minimum problem 

(1.4) mini/ | \7v\ 2 dx : v £ L^ 2 (Sl\K) , v = g on d D Sl\K\ . 

v L Jn\K > 

In these lectures, given a function u £ L 1 ' 2 (Q,\K) for some K £ /C(O), we always extend 
Vu to ft by setting Vu = a.e. on K. Note that, however, V« is the distributional 
gradient of u only in Q\K , and, in general, it does not coincide in Q with the gradient of 
an extension of u . 

To study the continuous dependence on K of the solutions of problem (1.2) we consider 
the Hausdorff distance between two sets K\, £ JC(fl) , which is defined by 

dn(Ki, K 2 ) := max { sup dist(x, K 2 ), sup dist(y, i^i)} , 

xeK! y£K 2 

with the conventions dist(x, 0) = diam(Q) and sup0 = 0, so that dn(&, K) = if K = 
and dn{0,K) = diam(fi) if K ^ 0. We say that (K n ) converges to K in the Hausdorff 
metric if dn(K n ,K) — > 0. The following compactness theorem is well-known (see, e.g., [30, 
Blaschke's Selection Theorem]). 

Theorem 1.1. Let (K n ) be a sequence in /C(O) . Then there exists a subsequence which 
converges in the Hausdorff metric to a set K £ /C(O) . 

The following example shows that the convergence of (K n ) to K in the Hausdorff metric 
does not imply the convergence of the solutions of problems (1.2) relative to K n to the 
solution relative to K , if we have no bound on the number of connected components of K n . 
In the next lecture we will prove a convergence result under a uniform bound on the number 
of connected components of K n . 

Example 1.2. Let Q. := (0, 1) x (-1, 1) , let d D n := (0, 1) x {-1,1} , and let g := ±1 on 
(0, 1) x {±1} . For every n > 1 let 

n — l . . 
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and let u n be the solution of problem (1.1) relative to K n and g. Then (K n ) converges 
in the Hausdorff metric to the set K := [0,1] x {0} and (u n ) converges in L 2 (Sl) to the 
solution u of the problem 



(1.5) 



Am = in SI, 



— — = on dj^Sl = dSl\doSl ■ 
ov 

u = g on duSl ■ 



Since this solution is given explicitly by u(x\,X2) ■= x 2 , we see that it does not satisfy the 
Neumann boundary condition du/du = on K . 

Proof. Let Sl ± := {x G SI : ±x 2 > 0}. By (1.3) we have / n± \Vu n \ 2 dx < J n \Vg\ 2 dx for 
every n, so that, passing to a subsequence, we may assume that (u n ) converges weakly in 
H 1 (Sl + l)Sl~) to a function w G H 1 (Sl + US1~) such that w = g on doSl. By symmetry 
the traces of u n from ST^ 1 vanish on K\K n . Let x„ be the characteristic function 
of K\K n . From the definition of K n it follows that (xn) converges to 1/2 weakly in 
L 2 (K,TL 1 ), where H 1 is the one-dimensional Hausdorff measure. Since the trace operator 
is compact from H 1 (Sl ± ) into L 2 (K, TL 1 ) , the traces of u n from SV^ converge to the 
corresponding traces of w strongly in L?{K,'H 1 }. Therefore (u^Xn) converges to w ± /2 
weakly in L X {K, TL 1 ) . As u^Xn = on if (recall that = on K\K n and Xn = 
on Kn), we conclude that w ± /2 = on K, therefore w G H 1 (S}). By using the weak 
formulation (1.2) we obtain 

/ Vu n Vzdx = Q VzeH 1 ^), z = 0ond D Sl. 
Passing to the limit as n — > oo we obtain 



/ 

Jn 



Vw\7zdx = Vz G H 1 ^) , z = on d D Sl . 

n 



This implies that w coincides with the solution of problem (1.5). By uniqueness, the whole 
sequence (u n ) converges to u. □ 

Remark 1.3. The hypothesis g := ±1 on [0, 1] x {±1} was intoduced only to simplify 
the proof. Indeed the same result holds when g is an arbitrary function of H 1 (Sl). To 
prove this fact, we can not use the equality u n = on K\K n , which is not true in the 
general case; instead we introduce the functions u n {x\,X2) ■= u n (xi,X2) — u n {xi,— x%) 
and w(x\,X2) := w{x\,X2) — w{x\, —X2), and observe that (u n ) converges to w weakly in 
H 1 (Sl + UO~). Since the traces of u n from Sl^ vanish on K\K n , arguing as before we 
obtain that the traces of w vanish on K . This implies that w + = w~ on K , and hence 
w G H 1 (S1). The conclusion follows now as as in the previous proof. 

In some cases the limit problem can contain a transmission condition, as shown by the 
following example, for which we refer to [16] and [27]. Note that in this case the one- 
dimensional measure of K n converges to the one dimensional measure of K . 

Example 1.4. Let SI := (0,1) x (-1,1), let d D Sl := (0,1) x {-1,1}, and let g G ff^fi). 
For every n > 1 let 

n—l . • 1 -I 

K n := |J [I,i±i- e - n l x{0}, 
w in n J 



i=0 



and let u n be the solution of problem (1.1) relative to K n . Then (K n ) converges in the 
Hausdorff metric to the set K := [0,1] x {0} and («„) converges in L 2 {Sl) to the solution 
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in Sl\K, 

on dfl\(dEtoUK), 
on doQ , 
u + ) on K , 

= {x e ft : ±X2 > 0} and v ± is the outer unit 

In the literature of homogcnization theory one can find other examples where the conver- 
gence in the Hausdorff metric of K n to K does not imply the convergence of the solutions 
of the Neumann problems on Q\K n to the solution of the Neumann problem on Q\K (see, 
e.g., [23], [1], and [14]). These papers show also that the bound on the number of con- 
nected components of K n , that we shall consider in the next lecture, would not be enough 
in dimension larger than two. 

The one-dimensional Hausdorff measure Ti 1 is not lower semicontinuous with respect to 
the convergence in the Hausdorff metric. For instance, in Example 1.2 we have T0-{K n ) = 
1/2 for every n, while for the limit set we have Ti}{K) = 1. However the lower semiconti- 
nuity holds if we have a uniform bound on the number of connected components, as shown 
by the following theorem. 

Theorem 1.5. Let (K n ) be a sequence in fC(Q) which converges to K in the Hausdorff 
metric. Assume that each set K n has at most m connected components. Then K has at 
most m connected components and 

V}(K C\U)< liminf H\K n n U) 

n — >oo 

for every open set U C M 2 . 

Proof. The case m = 1 and U = M 2 is the Golab theorem, for which we refer to [19, 
Theorem 3.18]. For an independent proof for an arbitrary U see [26, Theorem 10.19]. The 
case m > 1 is an easy consequence of the case m = 1 (see [15, Corollary 3.3]). □ 

2. CONVERGENCE OF SOLUTIONS 

Given an integer m > 1 and a constant A > 0, let /C^(O) be the set of all compact 
subsets K of f2, with "H (K) < A, having at most m connected components. 

In this lecture we give the main ideas of the proof of the following theorem. The complete 
proof can be found in [15, Section 5]. 

Theorem 2.1. Let m > 1 and A > 0, let (K n ) be a sequence in /C„(0) which converges 
to K in the Hausdorff metric, and let (g n ) be a sequence in i/ 1 (Sl) which converges to g 
strongly in H 1 ^). Let u n be a solution of the minimum problem 

(2.1) mini / \Vv\ 2 dx : v £ L 1>2 (0\K„) , v = g n on d D fl\K n \ , 

v 1 Jn\K n > 

and let u be a solution of the minimum problem 

(2.2) mini/ | Vw| 2 dx : v e L^^R) , v = g on d D fl\K) . 

v L JQ\K > 

Then Vu n — > Vu strongly in L 2 (Q;R 2 ) . 



u of the problem 



(1.6) 



Aw, = 

£ = 
ov 

u = g 
9u± 



dv ± 



where vr 1 is the restriction of u to il ± := 
normal to il ± . 
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This result is related to those obtained by A. Chambolle and F. Doveri in [12] and by 
D. Bucur and N. Varchon in [8], [9], and [10], which deal with the case of a pure Neumann 
boundary condition. Since we impose a Dirichlet boundary condition on do^\K n and a 
Neumann boundary condition on the rest of the boundary, our results can not be deduced 
easily from these papers, so we give an independent proof, which uses the duality argument 
which appears also in [10]. 

To focus on the main ideas of the proof, we consider only the case m = 1 and Cl simply 
connected. Moreover we assume that K PI dCl = , to avoid minor difficulties arising at 
the boundary. The technicalities needed to avoid these simplifying hypotheses can be found 



First of all, we want to construct the harmonic conjugate of u„ . Let R be the rotation 
on R 2 defined by , 2/2) := (—3/2, 2/1) ■ 

Definition 2.2. We say that a function v G i? 1 (SI) is equal to a constant c on a set 
K G K(Cl) if there exists a sequence (v n ) in C 1 (fi) converging to v strongly in if 1 (SI) and 
such that each v n is equal c in a neighbourhood of K . 

Remark 2.3. It is possible to prove that, if 



then v = c on K in the sense of the previous definition (it is enough to adapt the proof of 
[7, Thcorcmc IX. 17] or to apply [22, Theorem 4.5]). 

Theorem 2.4. Let K be a connected compact set contained in Cl and let u be a solution 
of problem (1.2). Then there exists a function v G i^ 1 (SI) such that Vv = RVu a.e. on 
Cl. Moreover v is constant on K and on each connected component of dj^Cl (according to 
Definition 2.2). 

Proof. If ipG C C °°(S1) , we have 

(2.3) VuVipdx = / VuVtpdx = 0, 



where the first equality follows from our convention Vu = a.e. in K , while the second one 
follows from (1.2). Equality (2.3) means that div(Vu) = in V'(Cl), hence rot(RVu) = 
in T>' (Cl) . As SI is simply connected and has a Lipschitz boundary, there exists v G H (Cl) 
such that Vu = R Vu a.e. in SI . 

Since du/dv = on OnCI, the tangential derivative of v (which is equal to the normal 
derivative of u) vanishes on S/vSl, and this implies that v is constant on each connected 
component of d^Cl. 

If K has a non-empty interior and a smooth boundary, then v is constant a.e. on 
the interior of K, since Vv = a.e. on K . Therefore v is constant on K according to 
Definition 2.2. 

The case of a general K G JC(d) can be obtained by approximating if by a decreasing 
sequence of compact sets with non-empty interior and a smooth boundary (we refer to [15, 
Theorem 4.2] for the details). □ 

Proof of Theorem 2.1. Note that u is a minimum point of (2.2) if and only if it satisfies 
(1.2); similarly, u n is a minimum point of (2.1) if and only if it satisfies (1.2) with K and 
g replaced by K n and g n . 

Taking v :— g n in the functional to be minimized, we obtain that the sequence (Vu n ) 
is bounded in L 2 (Cl;R 2 ). Therefore, passing to a subsequence, (Vu n ) converges weakly in 
L 2 (S1;R 2 ) to a function tp . As K G /C^(O) by Theorem 1.5, and hence meas(-ff) = 0, is 
easy to see that there exists a function u* G L 1 ' 2 (Cl\K) , with u* = g on <9d£1, such that 
tp = Vu* a.e. on Cl (we are assuming here that K n drjCl = 0, sec [15, Lemma 4.1] for the 
details). 



in [15]. 



lim ess sup \v — c\ = V.t G K , 

B T {x) 
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We will prove that 
(2.4) Vu* = Vu a.c. in Q\K . 

As the limit does not depend on the subsequence, this implies that the whole sequence 
(Vu„) converges to Vm weakly in L 2 (f2;R 2 ). Taking u n — g n and u — g as test functions 
in the equations satisfied by u n and u , we obtain 

/ \Vu n \ 2 dx = j Vu n Vg n dx, / \\Iu\ 2 dx = I SJuVgdx. 
Jn Jn Jn Jn 

As Vu n ^Vu weakly in L 2 (fl; M. 2 ) and Vg n — > Vg strongly in L 2 (Q; M. 2 ) , from the previous 
equalities we obtain that || Vu n \\L 2 (n;R 2 ) converges to || Vu\\ 1,1 (n-m?) , which implies the strong 
convergence of the gradients in L 2 (f2; M 2 ) . 

By the uniqueness of the gradients of the solutions, to prove (2.4) it is enough to show 
that u* is a solution of (1.2). This will be done by considering, for each u n , its harmonic 
conjugate v n given by Theorem 2.4. By adding a suitable constant, we may assume that 
J n v n dx = for every n. Since Vv n — R\7u n a.e. on ft, we deduce that (Vv n ) converges 
to RVu* weakly in L 2 (f2;R 2 ), and by the Poincare inequality (v n ) converges weakly in 
i^ 1 (SI) to a function v which satisfies Vv = RVu* a.e. on ft. 

Let us prove that v is constant on K according to Definition 2.2. This is trivial if K 
reduces to a point. If K has more than one point, then liminf„ diam(if„) > 0; since the sets 
K n are connected, we obtain also liminf„ cap(/4T„, ft) > 0, where the capacity cap(if„,£l) 
of K n with respect to is defined by 



cap(X„,fi) := inf j y \Vz\ 2 dx : z e C*(fi) 



z > 1 on K„ 



As v n = c n on K n for suitable constants c„ (see Theorem 2.4), using the Poincare inequality 
(see, e.g., [33, Corollary 4.5.3]) it follows that (v n — c n ) is bounded in i? 1 (fi), hence the 
sequence (c„) is bounded, and therefore, passing to a subsequence, we may assume that 
(c„) converges to a suitable constant c. 

Let us fix a constant R > with R < diam(if)/2 and R < dist(K,dfl) . Since Av n = 
on Sl\K n , by Maz'ya's estimate (see [24, Theorem 1]) there exist two constants M > and 
P > 0, independent of n, x n , and r, such that for every x n G and every r G (0, R) 

(2.5) esssup \v n ~ c n \ < Mexp ( - (3 / 7„(x„,p) — ), 

B r (a: n ) V Jr P ' 

where J n (xn,p) ■— ca,p(K n n B p (x n ), B 2 p(x n )) ■ For n large we have diam(if„) > 2_R, 
so that for every x n G X n and every p G (0, i?) we have i^n fl dB p (x n ) 7^ 0. As 
is connected, there exists a constant a > , independent of n, x n , and p , such that 
Jn(x n ,p) = cap(K n C\B (x n ),B 2 p{xn)) > a. Therefore (2.5) yields 

/ r \ a P 

ess sup \v n - c\ < M I — + |c„ - c\ 

for every x n G and every r G (0, i?) . 

Let us fix x G K and a sequence x n G K„ converging to x. For every p G (0, r) we have 
B p (x) C B r (x n ) for n large enough, hence 

/ r s a/3 

esssup \v n - c| < M — + |c„ - c| 
B p (i) V -K y 

for n large enough. Passing to the limit first as n — > 00 and then as p — > r we get 

(2.6) esssup \v - c\ < M (^-Y^ . 

B r (x) v -^ y 

As r — > we obtain that v is equal to c on K (sec Remark 2.3). 
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On the other hand, every v n is constant on each connected component of 8nQ (see 
Theorem 2.4). Since v n — weakly in we conclude that v is constant on each 

connected component of <9jvfL 

Therefore there exists a sequence (w n ) in C 1 (ri) converging to v strongly in H 1 ^), 
such that each w n is constant in a neighbourhood of K and in a neighbourhood of each 
connected component of 9/vfL 

Let z G L 1,2 (f2\A) with z = on do^-, and let ip n G C 1 (il) with tp n = 1 on supp(Vu>„) 
and (p n = on a neighbourhood of KuOnCI. As div(i? Vw„) = in £1 and z <p„ G #o(^) > 
we have 

/ RVw n Vzdx = / i? Vw n V(z(,5„) dx = . 
>/n\if Jn\K 

Since i? Vu* = Vw a.e. on , passing to the limit as n — > oo we obtain 



/ Vu*Wzdx = - RVvVzdx = 0, 
Jn\K Jn\K 



m\K Jn\K 
which shows that u* is a solution of (1.2). □ 

3. A QUASI-STATIC MODEL FOR BRITTLE FRACTURES 

Since the pioneering work of A. Griffith [21], the growth of a brittle fracture is considered 
to be the result of the competition between the energy spent to increase the crack and 
the corresponding release of bulk energy. This idea is the basis of the celebrated Griffith's 
criterion for crack growth (see, e.g., [31]), and is used to study the crack propagation along 
a preassigned path. The actual path followed by the crack is often determined by using 
different criteria (see, e.g., [18], [31], [32]). 

Recently G.A. Francfort and J. -J. Marigo [20] proposed a variational model for the quasi- 
static growth of brittle fractures, based on Griffith's theory, where the interplay between 
bulk and surface energy determines also the crack path. 

The purpose of this and of the next lecture is to give a precise mathematical formulation 
of a variant of this model in the two-dimensional case, and to prove an existence result for 
the quasi-static evolution of a fracture by using the time discretization method proposed 
in [20]. 

To simplify the mathematical description of the model, we consider only linearly elastic 
homogeneous isotropic materials, with Lame coefficients A and \x. We restrict our analysis 
to the case of an anti-plane shear, where the reference configuration is an infinite cylinder 
OxR, with fi C K 2 , and the displacement has the special form (0, 0, u(x\, X2)) for every 
(xi,X2,y) G fixR. We assume also that the cracks have the form KxR, where AT is a 
compact set in fl . In this case the notions of bulk energy and surface energy refer to a finite 
portion of the cylinder determined by two cross sections separated by a unit distance. The 
bulk energy is given by 

(3.1) £ / \Wu\ 2 dx, 

1 Jn\K 

while the surface energy is given by 

(3.2) k7i l {K), 

where k is a constant which depends on the toughness of the material, and H 1 is the one- 
dimensional Hausdorff measure, which coincides with the ordinary length in case A is a 
rectifiable arc. For simplicity we take fi = 2 and k = 1 in (3.1) and (3.2). 

As in the previous lectures is a fixed bounded connected open subset of E 2 with Lipschitz 
boundary. Following [20], we fix an open subset c^fi of <9f2, on which we want to prescribe a 
Dirichlet boundary condition for the displacement u. As in the previous lectures we assume 
that do^t has a finite number of connected components. 
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Given a function g on <9d^, we consider the boundary condition u = g on 8d^\K . We 
can not prescribe a Dirichlet boundary condition on Qofi n K , because the boundary dis- 
placement is not transmitted through the crack, if the crack touches the boundary. Assuming 
that the fracture is traction free (and, in particular, without friction), the displacement u 
in Q\K is obtained by minimizing (3.1) under the boundary condition u = g on djjfl\K . 
The total energy relative to the boundary displacement g and to the crack determined by 
K is therefore 

(3.3) £{g,K) = mini / \Vv\ 2 dx + H 1 (K) : v e L 1 ' 2 (ft), v = g on d D n\K) . 

v 1 Jq\k > 

The existence of a minimizer has been proved in the first lecture. 

In the theory developed in [20] a crack with finite surface energy is any compact subset 
K of H with H l (K) < +oo. For technical reasons, due to the hypotheses of Theorem 2.1, 
we propose a variant of this model, where we prescribe an a priori bound on the number of 
connected components of the cracks. Without this restriction, some convergence arguments 
used in the proof of our existence result are not justified by the present development of 
the mathematical theories related to this subject. Given an integer m > 1 , let /C^(£l) be 
the set of all compact subsets K of f2, with 11} (K) < +oo, having at most m connected 
components. 

We begin by describing a discrete-time model of quasi-static irreversible evolution of a 
fracture under the action of a time dependent boundary displacement g(t) , < t < 1. As 
usual, we assume that g(t) can be extended to a function, still denoted by g(t) , which 
belongs to the Sobolev space if 1 (ft) . For simplicity, we assume also that g(0) = 0. 

Given a time step 8 > 0, for every integer i > we set if := iS and gf := g(tf). 
The fracture Kf at time t\ is defined inductively in the following way: for i = we set 
Kq := Kq , while for i > 1 Kf is any minimizer of the problem 

(3.4) min {£(gf, K) : K G }C f m (Q), K D KU } . 
Lemma 3.1. There exists a solution of the minimum problem (3.4). 

Proof. By definition := Kq G JC^Q). A ssumc by induction that Kf_ 1 G JC^SY) and 
let A be a constant such that A > £{gf,Kf_ 1 ). Consider a minimizing sequence (K n ) 
of problem (3.4). We may assume that K n e /C^(Q) for every n. By the Compactness 
Theorem 1.1, passing to a subsequence, we may assume that (K n ) converges in the Hausdorff 
metric to some compact set K containing Kf_ 1 . For every n let u n be a solution of the 
minimum problem (3.3) which defines £(gf, K n ) . By Theorem 2.1 (Vw„) converges strongly 
in L 2 (Q;M. 2 ) to Vu, where u is a solution of the minimum problem (3.3) which defines 
£{gf,K). By Theorem 1.5 we have K G K m (ty and H^K) < liminf„ H 1 ^) < A, hence 
K G /C^(TT). As ||Vu|| = lim„||Vu„||, we conclude that £(gf,K) < limm£ n £(gf,K n ). 
Since (K n ) is a minimizing sequence, this proves that K is a solution of the minimum 
problem (3.4). □ 

Let uf be a solution of the minimum problem (3.3) which defines £(gf,Kf). Then the 
pair (uf, Kf) minimizes the sum of the bulk and surface energy among all K G /C/„(f2) with 
K D Kf_ x and among all u G L 1,2 (fl\K) with u = g on djjfl\K . 

In order to pass to a continuous-time model, given T > we define the step functions 
K S : [0,T] - /C4(0) and u s : [0, T] -> L 2 oc (Q) by setting 

(3.5) K s (t):=Kf and u s (t) := uf for t\ < t < tf +1 . 

Our purpose is to pass to the limit as S — > 0. To this aim we use the following result, which 
is the analogue of the Hclly theorem for compact valued increasing functions. 

Theorem 3.2. Let (K n ) be a sequence of increasing functions from [0,T] into K.(Q) , i.e., 
K n (s) C K n (t) for < s < t < T . Then there exist a subsequence, still denoted by (K n ) , 
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and an increasing function K: [0,T] — > /C(O), smc/i that K n (t) — > AT(i) m f/ie Hausdorff 
metric for every t G [0, T] . 

To prove Theorem 3.2 we use the following result, which extends another well known 
property of real valued monotone functions. 

Lemma 3.3. Let K\, K 2 : [0, T] — > /C(f2) 6e two increasing functions such that 

(3.6) ifi(s) C JT 2 (t) and # 2 (s) C 

/or every s, t G [0, T] witft s <t. Let be the set of points t G [0, T] such that K\(t) = 
K 2 (t) . Then [0, T]\0 is at most countable. 

Proof. For i = 1,2, consider the functions O x [0, T] — > K defined by f%(x,t) := 
dist(x, Ki(t j) , with the convention that dist(a;, 0) = diam(fi) . Then the functions fi(-,t) 
are Lipschitz continuous with constant 1 for every t G [0, T] , and the functions fi(x, •) are 
non-increasing for every x G ft . 

Let D be a countable dense subset of Q . For every x E D there exists a countable set 
N x G [0, T] such that fi(x, •) are continuous at every point of [0,T]\N x . By (3.6) we have 
fi(x,s) > f 2 (x,t) and f2(x,s) > fi(x,t) for every x G f2 and every s, t G [0,T] with s <t. 
This implies that fi(x, f) — f 2 (x, t) for every x G D and every t G [0, T]\N X . Let AT be the 
countable set defined by N := \J xeD N x , and let t G [0,T]\AT. Then ,/i(>,i) = f 2 {x,t) for 
every x € D , and, by continuity, for every x G f2 , which yields ATi (t) = if 2 (t) ■ This proves 
that [0, T]\N C 6, hence [0,T]\6cJV. □ 

Proo/ 0/ Theorem 3.2. Let D be a countable dense subset of (0, T) . Using the Compactness 
Theorem 1.1 and a diagonal argument, we find a subsequence, still denoted by (K n ) , and 
an increasing function K : D — > /C(f2) , such that K n {t) —> K(t) in the Hausdorff metric for 
every i G D . Let if" : (0,T] /C(H) and AT+ : [0,T) /C(O) be the increasing functions 
defined by 

Jf-(t):=cl( (J #(*)) for0<t<T, 

s<t, sec 

K+(t):= p| AT(s) for0<t<T, 
s>t, seD 

where cl denotes the closure. Let O be the set of points t G [0, T] such that K~(t) = K + (t) . 
As K~ and K + satisfy (3.6), by Lemma 3.3 the set [0,T]\O is at most countable. 

Since K~(t) G K(t) C K+(t) for every t G £», we have i<T(i) = AT"(t) = K+(t) for 
every f G 6 n D. For every t G 6\D we define AT(t) := K~{t) = K+{t). To prove that 
K n (t) — > AT(t) for a given t G 6\Z) , by the Compactness Theorem 1.1 we may assume 
that K n (t) converges in the Hausdorff metric to a set K* . For every s\, s 2 G D, with 
Si < t < s 2 , by monotonicity we have AT(si) G K* C AT(s 2 ). As AT* is closed, this implies 
AT"(t) G K* G K + (t) , therefore AT„(t) AT(t) by the definitions of 6 and K(t) . 

Since [0, T] \ (6 U D) is at most countable, by a diagonal argument we find a further 
subsequence, still denoted by (K n ), and a function K: [O,T]\(0Uf) — > /C(Q), such that 
AT„(t) -> AT(t) in the Hausdorff metric for every t G [0, 7]\(9UD). 

Therefore K n {t) — > AT(t) for every t G [0, T] , and this implies that AT(-) is increasing on 
[0,T]. □ 

The following result on the continuity of compact valued increasing maps will be used in 
the next lecture. Its proof is similar to the proof of Theorem 3.2. 

Proposition 3.4. Let K : [0,T] — > K,(Cl) be an increasing function, and let K~ : (0, T] — > 
/C(f2) aKc? AT + : [0,T) — > JC(Cl) be the functions defined by 

(3.7) K-(t) := cl(\J s<t K(s)) for0<t<T, 

(3.8) K+(t):=n s>t K(s) for0<t<T, 
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where cl denotes the closure. Then 

(3.9) K-(t)cK(t)cK+(t) forO<t<T. 

Let 9 be the set of points t G (0, T) such that K~(t) = K+(t). Then [0,T]\9 is at most 
countable, and K{t n ) — ► K{t) in the Hausdorf metric for every t 6 9 and every sequence 
(t n ) in [0,T] converging to t. 

Proof. It is clear that K + (-) and K~(-) are increasing and satisfy (3.6). Therefore [0, T]\Q 
is at most countable by Lemma 3.3. 

Let us fix t G 9 and a sequence (t n ) in [0, T] converging to t. By the Compactness 
Theorem 1.1 we may assume that K(t n ) converges in the Hausdorff metric to a set K* . 
For every s\, s 2 G [0, T], with si < t < s 2 , we have if(si) C K{t n ) C if(s 2 ) for n large 
enough, hence if(si) C K* C K(s 2 ). As is closed this implies K~(t) C K* C K + (t), 
therefore K* = K(t) by (3.9) and by the definition of 9. □ 

According to Theorem 3.2, there exist a sequence (6k) converging to and an increasing 
function K : [0, T] -> /C(fi) such that, for every t G [0, T] , K 5 (t) -> in the Hausdorff 

metric as J tends to along this sequence. In the next lecture we will prove the main 
properties of the function K : [0,T] — > /C(£l) obtained in this way, which represents the 
continuous-time evolution of the fracture. To simplify the notation, when we write 6 — > 
we always mean that 6 tends to along the sequence (6k) considered above. 

4. PROPERTIES OF THE CONTINUOUS-TIME MODEL 

In this final lecture we prove some properties of the function K : [0, T] — > /C(O) defined 
as the limit, as 6 — > , of the functions K$ : [0, T] — ► fC(Cl) given by the discrete-time model. 
This function represents the quasi-static evolution of the crack in our continuous-time model. 

To prove these results, we assume that the function t i— » g(t) , which gives the imposed 
boundary displacement on 8d^ 1 is absolutely continuous from [0,T] into if 1 (fi). Its time 
derivative is a Bochner integrable function from [0,T] into iJ 1 (fi), which will be denoted 
by g(t) . For the main properties of absolutely continuous functions with values in a Hilbcrt 
space we refer, e.g., to [6, Appendix]. 

We begin with a crucial estimate for the solutions of the discrete-time problems. Here 
and in the rest of the lecture (-|-) and || • || denote the scalar product and the norm in 
L 2 (ft;K 2 ). 

Lemma 4.1. There exists a positive function u>(6) , converging to zero as 5 — * 0, such that 

(4.1) ||V^|| 2 + W 1 (^)<||V U f|| 2 +W 1 (^) + 2 f J (Vu s (t)\\7g(t))dt + u;(5) 

Jtf 

for 0<i < j with < T. 

Proof. Let us fix an integer r with i < r < j . From the absolute continuity of g we have 

t 6 

9r+l-9 5 r= [ r+1 g(t)dt, 

where the integral is a Bochner integral for functions with values in iJ 1 (fi). This implies 
that 

(4.2) V 5r 5 +1 -V 5r 5 = ( r+ \g(t)dt, 

Jt 5 r 

where the integral is a Bochner integral for functions with values in L 2 (fi; R 2 ) . 

As u s r + g s r+1 -g s r £ L^(n\K S r ) and u s r + g s r+1 - g 5 r = g s r+1 on d D n\K S r , we have 

(4.3) £(g 5 r+1 , K s r ) < ||V«f + Vg s r+1 - V^|| 2 + H 1 ^.) . 
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< liv^ll 2 

where 



By the minimality of u b r+l and by (3.4) we have 

(4.4) || V?4 +1 || 2 + U\K 5 r+1 ) - £(g s r+1 ,K s r+1 ) < £(g s r+1 ,K s r ) . 
From (4.2), (4.3), and (4.4) we obtain 

WVu^f+HHK^) < ||V«* + V 5 * +1 - V^IP+W 1 ^ 4 ) < 

<\\v4f +n\K s r ) + 2 j^\vui\vm)dt + \\vm\\dt) 2 < 

t 6 t 6 
H 1 (K s r ) + 2 [ r+1 (Vus(t)\Vg(t))dt + a(5) [ " +1 || Vg(t)\\ dt , 

a(6):= max / " + ' || Vg(t) \\ dt — >0 

0<r, t£<T J t s 

by the absolute continuity of the integral. Iterating now this inequality for i < r < j we get 
(4.1) with lu{5) := a(6) J* \\Wg(t)\\ dt. □ 

Lemma 4.2. There exists a constant X, depending only on g and K , such that 

(4.5) ||Vuf||<A and H 1 ^) < X 
for every S > and for every i > with t\ < T . 

Proof. As v := gf is admissible for the problem (3.3) which defines £(gf,Kf), by the 
minimality of uf we have ||Vuf || < \\Vgf\\ , hence \\Vu s (t)\\ < ||V5*(t)|| for every t G [0, T] . 
As t i— » g(t) is absolutely continuous with values in 1 (SI) , the function t i— ► ||V<7(t)|| is 
integrable on [0,T] and there exists a constant C > such that ||Vg(i)|| < C for every 
t G [0, T]. This implies the former inequality in (4.5). The latter inequality follows now 
from Lemma 4.1 and from the inequality ||Vuq|| 2 + H}(Kq) < Ti}(K(f) , which is an obvious 
consequence of the minimality of Uq and of the fact that <?q = <?(0) = 0. □ 

Lemma 4.3. Let X be the constant in Lemma 4-2. Then K$(t) G /C^ (SI) and K(t) G 
IC^Tl) for every t G [0, T] . 

Proof. By Lemma 4.2 we have ft 1 < A for every t G [0,T] and every 6 > 0. By 

Theorem 1.5 this implies K(t) G (ft) for every t G [0, T] . □ 

For every t G [0, T] let u(i) be a solution of the minimum problem (3.3) which defines 
£(g(t),K(t)). 

Lemma 4.4. For every t G [0,T] we /mue Va^i) — > Vu(i) strongly in L 2 (S1;R 2 ). 

Proof. As M^(i) is a solution of the minimum problem (3.3) which defines £(g$(t), K$(t)) , 
and #,5 (t) — > <?(t) strongly in 7? 1 (S1), the conclusion follows from Theorem 2.1 and Lemma 
4.3. □ 

The following lemma shows the minimality of the set K(t) for the functional £(g(t), •) 
with respect to sets K containing K(t) . 

Lemma 4.5. For every t G [0, T] we have 

(4.6) £(g(t),K(t))<£(g(t),K) VK e £>(fi) , K D K(t) . 

Proof. Let us fix t G [0, T] and K G IC{ n (U) with K D K(t) . Since K s (t) converges to K(t) 
in the Hausdorff metric as 6 — > 0, it is possible to construct a a sequence (Kg) in /C^(S1) , 
converging to K in the Hausdorff metric, such that Kg D Kg(t) and Ji}(Kg \Kg(t)) — ► 
T-^iKXKit)) as <5 ->■ 0. By Lemma 4.2 this implies that W 1 ^) is bounded as <5 -> 0. 
The main difficulty in the construction of Kg is the constraint on the number of connected 
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components. The proof of the details is quite long, but elementary, and is given in [15, 
Lemma 3.5]. 

Let vg and v be solutions of the minimum problems (3.3) which define £(gg(t), Kg) 
and £(g(t),K), respectively. By Theorem 2.1 Vvg -> Vv strongly in L 2 (Q;R 2 ). The 
minimality of Kg(t) expressed by (3.4) gives £(gg(t), Kg(t)) < £(gg(t), Kg) , which implies 
||Vu,5(i)|| 2 < \\Vv s \\ 2 + H 1 (Kg\Kg(t)). Passing to the limit as 5 -> and using Lemma 4.4 
we get ||Vw(t)|| 2 < ||V-y|| 2 + U 1 {K\K{t)). Adding W 1 ^*)) to both sides we obtain 

(4.6) . □ 

We can now pass to the limit in Lemma 4.1. 
Lemma 4.6. For every s, t with < s < t < T 

(4.7) \\Vu(t)\\ 2 + n\K{t)) < \\Vu(s)\\ 2 + U\K{s)) + 2 / (V«(r)|Vfl(r))dr . 



Proof. Let us fix s, t with < s < t < T . Given S > let i and j be the integers such that 
t\ < s < tf +1 and tj < t < t S j +1 . Let us define sg := tf and tg := tj. Applying Lemma 4.1 



we obtain 

(4.8) \\Vug(t)\\ 2 +n 1 (Kg(t)\Kg( s )) < ||V U5 ( S )|| 2 + 2 / (Vua(T)|Vfl(r))dT + w(J), 



with u)(5) converging to zero as 5 — > 0. By Lemma 4.4 for every r 6 [0,T] we have 
Vug(r) — > Vu(t) strongly in L 2 (il,M 2 ) as 5 — > 0, and by Lemma 4.2 we have ||Viti(r)|| < A 
for every r G [0, T] . 

Given e > 0, let K £ (s) := {x e O : dist(x, AT(s)) < e}. As AT^s) C A" £ (s) for 5 small 
enough, we have Kg(t)\K e (s) C Kg(t)\Kg(s). Applying Theorem 1.5 with [/ = K 2 \A: £ (s) 
we get 

W 1 (A"(i)\A" £ (s)) < liminf ^(^(t^A^s)) < liminf H 1 (Kg(t)\Kg(s)) . 

s — >o s — >0 

Passing to the limit as e — > we obtain 

W 1 (A"(t)\A'( S )) < liminf H^Kgit^Kgis)) . 

s — >o 

Passing now to the limit in (4.8) as 5 — > we obtain (4.7). □ 

We are now in a position to prove the absolute continuity of the function t £(g(t), K(t)) 
and to compute its derivative. 

Lemma 4.7. The function 1 1— » £(g(t), K(t)) is absolutely continuous on [0,T] and 

(4.9) j f £(g(t), K(t)) = 2(V«(t)|Vfl(t)) for a.e. t e [0, T] . 
Proof. Let < s < t < T . From the previous lemma we get 

(4.10) £(g(t), K{t)) - £(g(s),K(s)) < 2 f (Wu(r)\Vg(r)) dr . 

J s 

On the other hand, by Lemma 4.5 we have £(g(s), K(s)) < £(g(s), K(t)) . It is easy to see 
that the Frechet differential d£(g,K) of £(g,K) (with respect to g) is given by 

(4.11) d£(g,K)h = 2 I \7u g Vhdx, 

Jn\K 

where u g is a solution of the minimum problem (3.3) which defines £{g,K). Therefore we 
have 

£(g(t),K(t))-£(g(s),K(t)) = 2 [ (Vu(r,t)|V$(r)) dr , 

J s 
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where u(r,t) is a solution of the minimum problem (3.3) which defines £(g(r), K(t)) . To- 
gether with the inequality £(g(s), K(s)) < £(g(s), K(t)) , this implies 

(4.12) £(g(t),K(t)) - £(g(s), K(s)) > 2 / (V«(t, *)|V$(t)) rfr . 



Since there exists a constant C such that ||Vu(r)|| < ||V<7(r)|| < C and || Vu(t, t) | < 
||Vg(r)|| < C for s < r < t, from (4.10) and (4.12) we obtain 

\£(g(t),K(t))-£(g(s),K(s))\ <2C [ ||V$(r)||dr, 



which proves that the function t i— ► £(g(t), K(t)) is absolutely continuous. 

As Vu(r,i) ->■ Vw(t) strongly in L 2 (fi;R 2 ) when r -> i, if we divide (4.10) and (4.12) 
by t — s, and take the limit as s — > i— we obtain (4.9). □ 

The result of the previous lemma can be expressed equivalently in the following way. 

Lemma 4.8. The function 1 1— » £(g(t), K(t)) is absolutely continuous on [0,T] and 

(4.13) ^£(g(t),K(s)) =0 for a.e. t G [0,T] . 

as s=t 

Proof. Let be the set defined in Proposition 3.4. By (4.11) and by Theorem 2.1 the 
differential d£(g,K(s)) tends to d£(g(t),K(t)) as (g,s) tends to (g(t),t) in iJ 1 ^) xl. 
Let us fix a point t in such that the function s i— > £(g(s), K(s)) is diffcrcntiablc at s = t 
and f is a Lebesgue point of g. For every s£ [0, T] we have 

= £(g(s),K(s)) - £(g(t), K(s)) + £(g(t), K(s)) - £(g(t),K(t)) = 
d£(g(T),K(s)) g(r) dr + £(g(t),K(s)) - £(g(t), K(t)) . 



Dividing by s — t and taking the limit as s — > t we obtain 



d_ 
ds' 



£(g(s),K(s)) = d£(g(t),K(t))g(t) + ±£(g(t),K(s)) 
-=t ds 



s=t 

The conclusion follows from (4.11) and Lemma 4.7. □ 

The properties of the function K: [0, T] — > /C(O) are summarized by the following theo- 
rem, which is an immediate consequence of Lemmas 4.3, 4.5, 4.7, and 4.8. 

Theorem 4.9. Let m > 1 , let g 6 4C([0, T]; H 1 ^)) , and let K G K f m (Q) . Then the 
function K: [0,T] — > /C(Q) introduced at the end of the last lecture satisfies the following 
properties: 

(a) A-(0)=tf , 

(b) K Q C if (s) C if (t) and K{t) G /C^(H) /or < s < f < T, 

(c) V0<t<T f(5(t),A-(t)) <£(g(t),K) V K e & m (Ti) , K D K (t) , 



(d) t i ► £(g(t),K(t)) is absolutely continuous on [0,T], 

(e) ^-£(g(t),K(s)) =0 /or a.e. t e [0, T] . 

as s=t 

Moreover every function K: [0,T] — > /C^(f2) which satisfies (a)-(e) satisfies also 

(f) j t £{g(t) 1 K{t)) = 2{Vu{t)\Vg(t)) for a.e. t G [0,T] , 
w/iere u(i) is a solution of the minimum problem (3.3) which defines £(g(t), K(t j) . 
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In our continuous-time model, the function K: [0,T] — ► tC^fl) represents the quasi- 
static irreversible evolution of the crack starting from Kq (condition (a)) under the action 
of the boundary displacement g(t). Condition (b) reflects the irreversibility of the evolution 
and the absence of a healing process. Condition (c) is a unilateral minimality condition. 
Condition (e) says that, for almost every t G [0, T], the total energy s i— ► £(g(t), K(s)) is 
stationary at s = t. Conditions (c) and (c) together lead to Griffith's analysis of the energy 
balance in our model, and, under some very mild regularity assumptions on the cracks, allow 
to express the classical Griffith's criterion for crack growth in terms of the stress intensity 
factors at the tips of the cracks (see [15, Section 8]). 

We underline that, although we can not exclude that the surface energy Ti}{K{t)) may 
present some jump discontinuities in time (sec [20, Section 4.3]), in our result the total energy 
is always an absolutely continuous function of time by condition (d) . 

If do^l is sufficiently smooth, we can integrate by parts the right hand side of (f) and, 
taking into account the Euler equation satisfied by u(t) , we obtain 

(4.14) ±£(g(t),K(t)) = 2 f ^g^dH 1 for a.e. t G [0, T] , 

at Jd D n\K(t) av 

where v is the outer unit normal to dCl. Since the right hand side of (4.14) is the power of 
the force exerted on the boundary to obtain the displacement g(t) on dn^\K(t), equality 
(4.14) expresses the conservation of energy in our quasi-static model, where all kinetic effects 
are neglected. 

The discrete-time model described in the previous lecture turns out to be a useful tool 
for the proof of the existence of a solution K(t) of the problem considered in Theorem 4.9, 
and provides also an effective way for the numerical approximation of this solution (see [5]), 
since many algorithms have been developed for the numerical solution of minimum problems 
of the form (3.4) (see, e.g., [2], [28], [29], [3], [11], [4]). 
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